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ABSTRACT 

We investigate how the typical dust extinction of Ha luminosity from a star-forming galaxy 
depends upon star formation rate (SFR), metallicity and stellar mass independently, using a 
sample of ^90,000 galaxies from Data Release 7 of the Sloan Digital Sky Survey (SDSS). 
We measure extinctions directly from the Balmer decrement of each source, and while higher 
values of extinction are associated with an increase in any of the three parameters, we demon- 
strate that the fundamental property that governs extinction is stellar mass. After this mass- 
dependent relationship is removed, there is very little systematic dependence of the residual 
extinctions with either SFR or metallicity, and no significant improvement is obtained from 
a more general parameterisation. In contrast to this, if either a SFR-dependent or metallicity- 
dependent extinction relationship is applied, the residual extinctions show significant trends 
that correlate with the other parameters. Using the SDSS data, we present a relationship to 
predict the median dust extinction of a sample of galaxies from its stellar mass, which has a 
scatter of ^ 0.3 mag. The relationship was calibrated for Ha emission, but can be more gen- 
erally applied to radiation emitted at other wavelengths. These results have important applica- 
tions for studies of high-redshift galaxies, where individual extinction measurements are hard 
to obtain but stellar mass estimates can be relatively easily estimated from long-wavelength 
data. 

Key words: dust, extinction, galaxies: evolution, galaxies: high-redshift, galaxies: ISM 
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1 INTRODUCTION 

Galaxies which are currently forming stars emit a large fraction of 
their luminosity at ultraviolet (UV) wavelengths, where the output 
from young, massive stars peaks. Some fraction of this radiation is 
absorbed and re-emitted through the Ha recombination line, and 
the intrinsic luminosity of both of these indicators scales directly 
withtherateof formation of massive stars (> 10 Mp,, with lifetimes 
less than ~ 20 Myr), giving a sensitive, well-calibrated p robe of the 
current star formation rate (SFR; e.g. lKennicutllll998[) . However, 
interstellar dust grains attenuate this short-wavelength radiation, re- 
ducing the observed flux and re-processing it into the far-infrared, 
so that only a small fraction of the intrinsic Ha or UV luminosity 
of a galaxy can be measured. Studies of star-forming galaxies are 
therefore hindered by uncertainties about the nature and amount of 
dust extinction which has occurred. 

The typical Ha extinction of a galaxy. Ana, can vary by 
several magnitudes and has been shown to d epend upon galaxy 
prop erties such as lumino s ity (e.g. [W ang & He ckmanI 199f ), SFR 
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(e.g. 'Hopkins etal.' 1 200 ll ; I Sullivan et ab 2001; 'Serta et al, 
Garn et al. 2009b), s tellar mass (e . g. [B rinchmann et alj 
Stasiiiska et al.l |2004| ; IPannella et afl l2009h and metallicity (e.g 
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lHeckmanetaLlll998l; ISoissier et aflbOO^ ; lAsari et alJlJOOTh . Ex- 
planations have been proposed for all of these effects: more mas- 
sive galaxies may have built up a larger dust reservoir with which 
to attenuate Ha radiation; more actively star-forming galaxies are 
likely to have larger and more dusty star-forming regions; more 
metal-rich galaxies have a greater dust-to-gas ratio, and therefore 
higher extinction. 

However, these studies have typically focused on one de- 
pendence at a time, and while they confirm that a single value 
of extinction should not be assumed for all galaxies, they do 
not address the question of whether stellar mass, SFR or metal- 
licity is more fundamental. Further complications are introduced 
through the known correlations between galaxy pr operties - more 
massive galaxies tend to be m ore metal-rich (e.g. iTremonti et al.l 
I2OO4I ; iKewlev & EUisonl 120081) and, amongst star- forming galax- 
ies, more massive galax ies also hav e a greater typical SFR (e.g. 
iBrinchmann et al.1 |2004|PI Recently iMannucci et alj ilQldt have 
proposed the existence of a 'fundamental metallicity relation', 
whereby galaxies define a tight surface in the 3-dimensional space 
of mass, SFR and metallicity. They also argue that this surface does 
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The mass-SFR congelation break s down for masses a bove those at which 
'SFR-downsizing' has set in (e.g. ICowie et alj|l996l) . whereby the most 
massive galaxies are no longer star-forming in the nearby Universe. 
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not evolve out to z ~ 2.5, and that previous observations that the 
trends between mass, metallicity and SFR evolve with red shift (e.g. 
lErb et alJllOO^ : iDaddi etlj]|2003l : iMaiolino et al.ll2008l) arise be- 
cause the high-redshift observations sample different regions of the 
surface. What is certainly clear is that the large number of compet- 
ing effects make it difficult to determine what drives the observed 
systematic variation in dust extinction between different galaxies. 

The most direct method of estimating extinction makes use of 
the strong wavelength-dependence of dust attenuation - by com- 
paring the observed flux ratio of the Ha and Hj6 lines to the ratio 
that is expected in the absence of dust (the 'Balmer decrement'), the 
amount of extinction at a given wavelength can be calculated. How- 
ever, this method is only practical at low redshift, given the relative 
weakness of the Hj3 line, and the fact that the Ha line is redshifted 
out of the optical band at z > 0.5. ff a suitable parameterisation for 
extinction could be calculated from galaxies in the local Universe, 
it could usefully be applied to the high-redshift population, where 
extinction cannot be estimated as easily. 

In this work we aim to disentangle the effects that stellar mass, 
SFR and metallicity have upon dust extinction, determine which of 
these parameters is the most fundamental, and define a relationship 
to estimate the typical extinction of a galaxy when measurements 
of the Balmer decrement are not available. In Section[2]we present 
the data, describe our methods for selecting a clean sample of star- 
forming galaxies, and describe the estimates of SFR, stellar mass, 
metallicity and extinction that we use. Section [3] presents the ob- 
served variation in Ha extinction with each combination of galaxy 
parameters, and our models that describe how to predict dust ex- 
tinction. In Section |4] we discuss these results, and their implica- 
tions. 

Throughout this work we assume a concordance cosmology 
of O.M = 0.3, fiA = 0.7 and Hq = 70 km s"' Mpc"^ 



2 DATA SAMPLE 

2.1 Selection of a star-forming sample 

In this work we make use of the superb spectroscopic informa- 
tion available from the Sloan Digita l Sky Survey, Data Release 7 
(SDSS DR7; lAbazaiian etaLll2009h . This dataset has been anal- 
ysed by a group from the Max Planck Institute for Astrophysics, 
and Johns Hopkins University (hereafter the MPA/JHU groupj^ - 
descriptions of the a nalysis pipeline for a prev ious SDSS data re- 
lease can be found in lBrinchmann et alj ( 12004) and lTremonti et al.i 
( 120041) . Spectroscopic redshifts and line fluxes (from a Gaussian fit 
to continuum-subtracted data, corrected for Galactic reddening) are 
provided for 927,552 sources at z < 0.7. We exclude from our anal- 
ysis the sources spectroscopically classified as a 'QSO', sources 
outside the redshift range 0.04 < z < 0.2, within the redshift range 
0.146 < z < 0.148 (see Section [2.2. It . or sources with an uncer- 
tain estimate of their redshift (Z_WARNING > 0), and also remove 
multiple entries of a number of objects from the dataset (where we 
define a multiple entry to consist of two catalogued sources with 
angular separation less than 5 arcsec, and Az < 0.001). 617,822 
sources are retained after these criteria have been applied. 

In this work we require a sample of star-forming galax- 
ies, from which we can measure an extinction and metal- 
licity from emission-line ratios. In order to distinguish be- 
tween star-forming galaxies and AGN, we use the 'BPT' 
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Figure 1. Lower panel: the number of sources above a given Ha S/N ratio. 
Upper panel: the fraction of these sources that have a detection of the H/3, 
[O III] or [N II] lines above a given S/N ratio, and the fraction with detec- 
tions of all three lines, with S/N of > 3 for the Hj3 and [N II] lines, and S/N 
of > 2 for the [O III] line. 
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( IBaldwin. Phillips & TerlevichI Il98lh emission line diagnostic, 
which compares the [O IIl]5007/H/3 and [N Il]6584/Ha line flux 
ratios. However, in order to use this diagnostic, all four emission 
lines must be detected with sufficient significance that accurate line 
flux ratios can be calculated. Note that the formal line flux uncer- 
tainties quoted in the IVIPA/JHU DR7 catalogue significantly un- 
derestimate the true errors ~ we have increased these uncertainties 
by the recommended factors of 2.473, 1.882, 1.566 and 2.039 for 
the Ha, Hj3, [O ill] and [N II] lines respectively, which were calcu- 
lated by the MPA/JHU group through comparisons of the derived 
line fluxes of sources that were observed multiple times. 

We perform a primary selection on the signal-to-noise (S/N) of 
the Ha line (we investigate possible biases of this selection method 
in Section 3.5), and measure the fraction of these sources that also 
have a detection of the other lines with a S/N of > 3. Fig.[T]shows 
how the number of sources above a given Ha S/N ratio varies, 
along with the fraction of these sources with detections of the other 
emission lines above a given significance. It is clear that the limiting 
factor in this analysis is obtaining a secure detection of the [O III] 
line - any sample with an Ha S/N threshold > 10 is essentially an 
[O IIl]-selected sample. 

For this analysis we set an Ha S/N requirement of 20, which 
includes 143,472 sources. We set an Hj3 S/N requirement of 3, 
which includes 99.97 per cent of these (i.e. we are not biasing the 
sample through our Hj8 cut). Of these sources, 99.86 per cent had 
a detection of the [N II] line with S/N > 3. In order to not bias 
our sample against weak [O III] emitters (which have low metal- 
licity; see Section |2^ . we reduce the [Olll] S/N threshold to 2. 
This cut includes 139,586 sources (97 per cent), in comparison 
to the 89 per cent which would h ave been selected above a S/N 
of 3. The .Kauffmann et al.l ( l2003br) criteria was then used to iden- 
tify sources on the BPT diagram which are definitely star-forming 
(120,650 galaxies; 86 per cent of the sample). 



2.2 Sample properties 

2.2.1 Extinction 

For each galaxy we measure the Ha to H/3 line flux ratio, 
Sna/S^a, and estimate an extinction at a particular wavelength A, 
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^A' using 



A, = 






lo: 



'SlO 



2.86 



Sua /S] 



H/J 



(1) 



where 2.86 is the intrinsic Ha/H/3 line flux ratio that is appropri- 
ate for Case B recombination, a tem perature of T = 10" ^ K and an 
electron density of iig = 10 cm^ l IBrocklehurstll 197 ih . We con- 
sider the accuracy of using the single value of 2.86 for all galaxies 
in Section [23l 

We use the lCalzetti et alj J200Cft dust attenuation law to cal- 
culate the values of ^;l = Ai/E{B — V) at the wavelengths of the 
Ha and Hj3 emission lines. This relationship was measured from 
observations of z ~ starburst galaxies, and has been shown to be 
appropriate for galaxies at redshifts up to --^ 0.8 (Gam et al. 2009b|). 
Galaxies which had an observed value of Sua /Sun < 2.86 were as- 
signed an extinction of mag (0.6 per cent of the sample). 

An uncertainty in the extinction for each galaxy was calcu- 
lated, using the Ha and H/3 flux uncertainties and standard error 
propagation methods; the median uncertainty is 0.2 mag. The ex- 
tinction errors were not found to correlate with any of the other 
properties of the sample, although a significant increase in the typ- 
ical error was found for galaxies within the redshift range 0.146 < 
z < 0.148 (median uncertainty of 0.34). This appears to be due to 
contamination of the redshifted H/3 line with the [O l]5577 sky 
emission line JOsterbrock & Martellll992h - we reject from our 
analysis any sources within this redshift range. 



2.2.2 Stellar mass 

Stellar mass estimates, M*, have been calculated for the DR7 
data by the MPA/JHU tearrfj, based upon fits to the total 
(CMODEL) SPSS phot omet ry and following t he philosophy of 
iKauffmann eTaP J2()03ah and lSalimetd] ( l2007h . We adopt these 
mass estimates, and we exclude from our analysis all galaxies 
(7 per cent of the remaining sample) where this estimate was not 
available. Note that the excluded galaxies have comparable proper- 
ties to those retained, so this exclusion should not bias the results 
in any way. 

Additional stellar mass estimates have been calculated from 
photometry measured within the 3 arcsec diameter SDSS fibres. 
A comparison of these two estimates allows us to define an effec- 
tive aperture correction factor, which corresponds to the fraction of 
the stellar mass (and hence broadly, the stellar light) that is located 
within the fibre - this correction factor can be used to identify and 
reject highly extended galaxies, where the fibre-based spectroscopy 
will not be representative of the properties of the who le galaxy. This 
is effe ctively the same correction as is discussed bv iHopkins et al.l 
1 120031) , who use the ratio of total and fibre-based r-band magni- 
tude s to correct the detec ted Ha emission of SDSS galaxies. 

iKewlev et al.l J2005l) discuss the potential effects of estimat- 
ing Balmer extinctions, metallicities and SFRs from spectroscopic 
observations made within a limited aperture, and conclude that an 
aperture covering factor of 20 per cent is sufficient to minimize sys- 
tematic differences between these quantities and the global values. 
They recommend a minimum redshift of z = 0.04 in order to obtain 
an average fibre-covering factor of greater than 20 per cent for a 
typical galaxy, which we have adopted for this analysis. We further 
reject galaxies at higher redshift that have a fibre-coverage factor 
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< 20 per cent, along with a small number of galaxies where a con- 
version between point spread function magnitudes and fibre mag- 
nitudes performed by the MPA/JHU group has been unsuccessful 
(J. Brinchmann, private communication), leaving 95,240 galaxies 
in the sample (85 per cent of the galaxies which had stellar mass 
information available). 



2.2.3 Star formation rate 

In order to estimate SFRs, we calculate the intrinsic Ha lumi- 
nosity for each galaxy, using the aperture and extinction-corrected 
Ha line flux. From this luminosity we calculate a SFR using the 
relation s hip in .Kennicutt (199 8.) . modified slightly to assume a 
iKroupal ( 120011) Initial Mass Function (IMF) that extends between 
0. 1 and 100 Mq (to agree with the IMF used in the MPA/JHU stel- 
lar mass calculations). The choice of IMF is unimportant in this 
analysis as Ha emission is produced mainly by high-mass stars, 
where most IMFs agree. Implicit in our calculation of SFRs from 
an aperture-corrected Ha flux is the assumption that the radial dis- 
tribu tions of Ha flux and stellar continuum are comparable - both 
Brinchmann et al.l J2004h and lKewley et al. (2005) find that this is 
an acceptable approximation, if at least 20 per cent of the light of a 
galaxy is contained within the fibre. 



2.2.4 Metallicity 

We use the 03N2 indicator jPettini & Pagell 20040 as our proxy for 
'metallicity' (denoting the gas-phase abundance of oxygen relative 
to hydrogen), where 



03N2.1og,ofpHp^V 
V [Nnl6584/Ha J 



(2) 



We use this indicator for three reasons: (i) due to the similarity in 
wavelengths for the Ha and [N II] lines, and the Hj8 and [O III] 
lines, this indicator is essentially independent of the effects of dust 
attenuation; (ii) the 03 N2 indicator h as a unique metallicity for 
each line flux ratio (see lPettini & Pagel 2004. for further details); 
(iii) the emission lines that are required to construct the 03N2 diag- 
nostic are the same as we have used to select star-forming galaxies 
from the BPT diagram, and we can therefore measure a metallic- 
ity for our entire sample, without introducing additional selection 
criteria. 

The conv e rsion to an (O/H) abundance comes from 
IPettini & Pagel ( l2004f) . 12 + logio(0/H) = 8.73 - 0.32 x 03N2, 
where they caution that the relationship may break down for low- 
metallicity galaxies with 03N2 > 1 . 9, or 1 2 + log i q (O/H) < 8 . 1 - 
only 0.2 per cent of our sample have these values of 03N2, and we 
reject these from our analysis. All (O/H) abundances can be con- 
verted to solar metallicities, using a value of I2 + logjQ(0/H)::i = 
8.66 ( Asplund et al. 2004), and the range 8.1 < 12 + logio(0/H) < 
9 corresponds to 0.28 < Z/Zr,, < 2.19. 



2.2.5 Summary 

The details of the selection cuts applied to the sample are sum- 
marised in Table[T] The vast majority of the selection cuts are either 
bias-free (e.g. cuts in redshift), or remove so few galaxies that any 
bias must be insignificant. The possible exceptions to this are the 
selection of star-forming galaxies using a signal-to-noise thresh- 
old in Ha, the removal of galaxies without stellar mass estimates, 
and the fibre-coverage cut. As noted in Section 2.2.2, the properties 
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Table 1. Details of the selection cuts applied to the SDSS sample. 



0.2 



Reason 


tigais retained % r 


emoved by cut 


Initial Sample 


927,552 


- 


Restrict to 0.04 < z < 0.20 


666,978 


28.09% 


Reject 0. 146 <z< 0.148 


659,737 


1.09% 


Reject uncertain z 


658,464 


0.19% 


Reject QSOs 


657,321 


0.17% 


Remove multiple entries 


617,822 


6.01% 


Require S/N (Ha)>20 


143,472 


76.78% 


Require S/N (Hj3)>3 


143,433 


0.03% 


Require S/N ([N1I])>3 


143,237 


0.14% 


Require S/N ([01I1])>2 


139,586 


2.55% 


Select as SF gals using BPT 


120,650 


13.57% 


Reject if no mass estimate 


112,008 


7.16% 


Reject if <20% fibre coverage 


95,240 


14.97% 


Restrict to 03N2 > 1.9 


95,041 


0.21% 



(ie. magnitudes, redshifts. Ha fluxes, emission line ratios, etc) of 
the galaxies excluded due to an absence of a mass estimate show 
a comparable distribution to those of the retained galaxies, so this 
cut is not expected to introduce any significant bias. The fibre cov- 
erage cut removes galaxies with larger angular sizes; these tend to 
be at the lower end of the redshift range, and to be amongst the 
more massive galaxies at those redshifts. This population of galax- 
ies is well-represented at higher redshifts in the sample (where their 
angular sizes are smaller, but they remain luminous enough to be 
detected), so again their exclusion is expected to introduce only 
minimal bias. We conclude that the only selection cut which has the 
potential to introduce a significant bias is that of the Ha signal-to- 
noise threshold, and we will investigate this further in Section 3.5. 
The final number of galaxies in our sample is 95,041, and 
Fig. |2] shows the distributions of M,^, SFR, metallicity, redshift, 
extinction and extinction uncertainty. The median stellar mass 
is lo'"'*" Mq, with an inter-quartile range of 10^**^ Mq, to 
jq1().47 f^^^^ jjjg median SFR is 6.1 Mg, yr"', with an inter- 
quartile range of 2.9 to 12.9 Mq yr^'. The median metallicity is 
12 + logio(0/H) = 8.72, or Z = 1.1 Zq; the metallicity distribu- 
tion is asymmetric with an inter-quartile range of 8.61 to 8.78, and 
a longer tail towards low metallicities. The sample was restricted to 
a redshift range of 0.04 < z < 0.2 and is mainly at z < 0. 1 (median 
= 0.085; inter-quartile range 0.067 to 0.113). The galaxies in the 
higher redshift tail predominantly have high masses and high SFR, 
but their extinction properties are comparable to those of lower red- 
shift galaxies of the same mass and SFR - there is no evidence 
for redshift evolution within the sample. The median Ha dust ex- 
tinction is 1.03 mag, close to the canonical value of 1 mag often 
assumed for Ha-selected galaxies. The median uncertainty in ex- 
tinction estimates is 0.2 mag. 



2.3 Accuracy of the extinction estimates 

A previous SDSS data release has been analysed by the MPA/JHU 
gro up in the same manner as the DR7 data which we use in this pa- 
per. ISrinchmann et al.l ( l2004f) have used these data to test the two 
assumptions that we have made, namely that the intrinsic Ha/Hj3 
ratio is equal to 2.86, and that the conversion between intrinsic 
Ha luminosity and SFR that we use is appropriate for all galax- 
ies. They conclude that both assumptions are suitable when used 
in combination with each other, although each is slightly incorrect 
- the Case B ratio depends upon temperature and density, varying 
between extremes of 2.72 (for T = 2x 10 K, rig = 10 cm^ ) and 
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Figure 2. The fractional distribution of (a) stellar mass; (b) star-formation 
rate; (c) metallicity; (d) redshift; (e) Ha extinction and (f) uncertainty in 
Ha extinction for our sample. 



3.03 (r = 5 X 10^ K, ;ie = 10^ cm"^), while the most metal-rich 
galaxies pr oduce a slight ly lower Ha luminosity for the same SFR. 
Brinchmann et al.l 1 12004 1) conclude that the combination of the two 
assumptions approximately cancel out when calculating the SFR, 
but there is the potential for a bias in our extinction estimates if 
the wrong intrinsic line ratio has been assumed. However, even the 
maximum error in our choice of Case B ratio would only lead to an 
error in extinction of ±0. 16 mag, which is smaller than the typical 
uncertainty in the individual estimates of extinction. 

In order to look for systematic offsets in our estimate of ex- 
tinction with other properties of a galaxy, we match our SDSS 
sample with the 1.4-GHz Faint Images of the Radio Sky at Twenty- 
cm (FIRS T;lBecker et alj 19951) survey, using the methodology de- 
scribed in iBest et al.l ( 120051) . 2404 galaxies in our sample have a 
unique counterpart in FIRST, and for each galaxy we calculate 
a radio-based SFR using the relationship in Bell (2003), assum- 
ing that the radio spectrum of each galaxy can be modelled by 
Sv o= V"°^ (see e.g. lGarn et al.l2008t) . This SFR indicator has been 
shown to agree with infrare d-based values for star-forming galax- 
ies at redshifts up to z = 2 JGarn et al.ll2009 a-). and gives an inde- 
pendent estimate of the intrinsic SFR, without being biased by the 
effects of dust extinction. 

The resolution of FIRST is 5 arcsec, so the radio flux will be 
sampling a slightly larger, but comparable region of each galaxy 
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Figure 3. The difference in extinction estimates between those calculated from the Balmer decrement, ABaimci-: and those from a comparison with radio-based 
SFRs, Arajio, as a function of (a) stellar mass; (b) metallicity and (c) redshift. Individual galaxies are denoted by dots, and binned data (median, ± Iff) by 
points with error bars. 



to the SDSS fibres. We therefore compare the radio-based SFRs 
to SFRs calculated from the observed Ha flux, without correcting 
for extinction or aperture size, in order to obtain an estimate of the 
Ha dust extinction for each galaxy, Ajadlo^ that is independent of 
the value calculated from the Balmer decrement, ABaimer- Fig.[3ja) 
shows the agreement between these two extinction estimates - data 
have been binned by stellar mass, and for each bin we calculate 
the distribution of the value of AA = Agaii^er — Afajjo. The median 
value of AA is taken to be our indicator of the 'typical' value, and 
the values of AA at the 16th and 84th percentiles (±1ct) to indicate 
the range that each bin covers. 

A systematic offset is seen for all bins, which is equivalent 
to having the radio SFRs systematically 30 per cent greater than 
the extinction-corrected Ha SFRs. This offset is likely to be due to 
the fact that the FIRST and SDSS observations are estimating their 
SFRs from different portions of a galaxy, although there may also 
be some contribution from a systematic uncertainty in the two SFR 
calibrations. However, no significant mass-dependent variation is 
seen between the two estimates of extinction. Fig.jSjb) and (c) show 
the equivalent plots for galaxies binned by metallicity and redshift 
- again, no significant variation is seen in the two extinction esti- 
mates. We do not show the equivalent plot for data binned by SFR, 
as this is affected significantly by sample bias due to the FIRST 
completeness limit of 0.75 mjy corresponding to a SFR limit that 
is significantly higher than that of the Ha selection. 

The results of this section lead us to conclude that our extinc- 
tion estimates are robust, and that our choice of a single value of 
the intrinsic Ha/Hj3 ratio should not introduce mass, metallicity or 
redshift-dependent trends into our analysis. 



is overlaid (median, ± la). A clear increase is seen in the median 
extinction of galaxies with a higher SFR, in a greement with previ- 
ous studies of the low-redshift Universe (e.g. iHopkins et alJuOOll ; 
ISullivanetal. 2001 : Berta et al. 2003) . 



iGam et al.l d2009bl) found that there was a well-defined rela- 
tionship between dust extinction and SFR for a sample of Ha- 
emitting galaxies located at z = 0.84. Note that Garn et al. (2009b|) 
calculated Ha extinction using an independent method (a compar- 
ison of Ha and 24-flm SFR indicators, rather than from observa- 
tions of the Balmer decrement). We overlay the lGam et alj ( l2009bl) 
relationship on Fig. |4la) - the relationship is in excellent agree- 
ment, with both the normalisation and the slope (despite hints of a 
slight change in the latter) consistent to within la. This confirms 
the conclusions of that work, that the dependence of dust extinction 
on SFR does not change significantly between z~ and z = 0.84 
(at least over the SFR range that was covered, 3 - 100 Mq yr^'). 
This result also confirms that there should be no significant dif- 
ference in the typical extinction relations for galaxies at different 
redshifts within our sample (0.04 < z< 0.2). 



We repeat this analysis to determine the variation in typical 
extinction with galaxy metallicity, which i s shown i n Fig. |4jb). 
As has been found i n previous studies (e.g. iHeckman et aljil998l : 
iBoissier et"ai] |2004 lAsari et alj l2007h . the typical extinction of 
a galaxy increases with metallicity and is approximately 1 mag 
for solar- metallicity galaxies, which have 12 + log[g(0/H) = 8.66 
JAsplund et alj2004h . 



3 RESULTS AND ANALYSIS 

3.1 Variation in Ha extinction with star formation rate, 
metallicity and stellar mass 

Fig. Eta) shows the relationship between SFR and extinction for 
the galaxies in our sample - the grey-scale illustrates the region 
of extinction / SFR space that is populated with galaxies. In or- 
der to quantify the observed trend, we separate our sample into 20 
logarithmically-spaced bins, selected by SFR. Within each bin, we 
calculate the distribution of measured Ha extinction values, which 



Fig. |4jc) shows the equivalent results, binning the sam- 
ple by stellar mass. The typical extinction of a galaxy in- 
crea ses with stellar mass, as has been found in previous studies 
(e.g. iBrinchmann et alj|2004l : IStasihska erZll2004l : IPannella et all 
2009), and rises steadily between 10^ and 10^' Mq. With the ex- 
ception of the outer bins of data (which contain fewer galaxies) the 
typical lC7 width of the distribution is --^0.3 mag, comparable to 
the distribution widths seen on the extinction / metallicity plot, and 
substantially smaller than the ~ 0.5 mag distribution width of the 
extinction / SFR plot. 
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Figure 4. The variation in Ha extinction with (a) star-formation rate; (b) 
metallicity; (c) stellar mass. The number of individual galaxies within each 
pixel is indicated by the grey-scale, and the binned data (median, ± 1 (T) by 
the solid and dashed lines respectively. Polynomial fits to the median are 
shown by the dotted lines - see Section [T4l for more details. For compar- 
ison, the best- fitting relat i onship between SFR and extinction for z = 0.84 
galaxies from lGam et al.l (2009bj) is shown on sub-figure (a), by the dash- 
dot and dotted lines (median and uncertainty respectively). 



3.2 Disentangling the relationships between Ha extinction 
and other galaxy properties 

We have shown that the typical extinction of a galaxy correlates 
with its SFR, metallicity and stellar mass. In order to disentangle 
the complex relationships between these properties, we separate 
our sample by one of these properties into six sub-samples, and 
bin each sub-sample by both of the remaining two variables. For 
each bin, we calculate the distribution of Ha extinction (median 
±1(7). The size of the SDSS DR7 dataset means that we can carry 
out this binning on a very fine scale - in the following analysis, any 
bins containing fewer than 20 galaxies are not shown for clarity, 
and contours illustrate the regions where bins contain more than 
100 or 1000 galaxies. 

We first take six sub-samples of galaxies with increasing 



Table 2. The covariance matrix for the star-forming galaxies in our sam- 
ple. All parameters have been normalised to have a distribution with zero 
mean and unit standard deviation. 'SFR' = logio(SFR / Mr, yr"'), 'Mass' 
= logio(M«/M0), 'Metallicity' = 12 + log,o(0/H), 'Extinction' = Ana- 





'SFR' 


'Mass' 


'Metallicity' 


'Extinction' 


'SFR' 


1.000 


0.734 


0.316 


0.612 


'Mass' 


0.734 


1.000 


0.717 


0.719 


'Metallicity' 


0.316 


0.717 


1.000 


0.572 


'Extinction' 


0.612 


0.719 


0.572 


1.000 



metallicity, and calculate how extinction varies with mass and SFR 
for each of the sub-samples. If changes in metallicity are the dom- 
inant factor responsible for variations in the extinction of differ- 
ent galaxies, then we would expect to see an approximately con- 
stant extinction in each sub-sample, for any mass or SFR. How- 
ever, Fig. [5] shows the results: while different metallicity galaxies 
populate different regions of the diagram, and the typical extinction 
of a galaxy does increase for the higher-metallicity sub-samples, 
there remains a significant variation in extinction within each plot 
(in the sense that extinction increases for higher masses or SFRs) 
that can not be explained purely through a metallicity dependence. 
The well-known correlation between stellar mass and SFR of star- 
forming galaxies can be seen in Fig.[5]through the direction of the 
contoured regions, and appears to be equally strong for galaxies of 
all metallicities. 

Fig.|6]shows the equivalent results, for sub-samples of the data 
selected by SFR, and extinction as a function of mass and metallic- 
ity. The tight correlation between mass and metallicity means that 
the width of the region that is populate d in each SFR sub-sample is 
relatively small. i Tremonti et al.l ( 120041) have identified that the cor- 
relation begins to break down at stellar masses above lO'"'^ Mq-,, 
with more massive galaxies not continuing to be more metal-rich 
- this effect can be seen clearly in Fig. [6] As for Fig. |5] a varia- 
tion in the extinction of a galaxy within a given sub-sample can be 
seen, implying that variations in SFR are not the dominant factor 
responsible for differences in extinction between galaxies. 

Fig. [7] shows the extinction as a function of SFR and metallic- 
ity, for sub-samples selected by their stellar mass. The same effect 
as before is apparent - different regions in the SFR / metallicity 
plane are populated by galaxies in each of the stellar mass sub- 
samples - and a weak anti-correlation between SFR and metallic- 
ity can be seen, for galaxies of a particular stellar mass. However, 
in contrast to the previous two figures, there is little variation in 
the median extinction of a galaxy at any populated point in each of 
the sub-samples. This implies that the majority of the variation in 
extinction between galaxies is driven by variations in stellar mass, 
rather than SFR or metallicity. 



3.3 Principal component analysis 

In order to quantify the relative importance of each of the pa- 
rameters that we are considering, we normalise each to have zero 
mean and unit standard deviation, and calculate the covariance ma- 
trix of the resulting variables. Each entry in this matrix measures 
the strength of the linear dependence between two variables, and 
ranges between +1 (a perfect correlation between the variables), 
(no correlation) and —1 (perfect anti-correlation). Table |2] reports 
the results: mass is strongly correlated with all the other variables, 
and SFR and metallicity are also strongly correlated with extinc- 
tion. 
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Figure 5. The variation in Ha extinction with stellar mass and star-formation rate, for sub-samples of galaxies having metallicity between 8.2 < 12 + 
logio(0/H) < 8.9. The grey-scale represents the median extinction within each bin, and contours illustrate the regions where bins contain more than 100 and 
1000 galaxies respectively. The metallicity range and number of galaxies in each sub-sample are listed above each figure. 
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Figure 6. The variation in Ha extinction with stellar mass and metallicity, for sub-samples of galaxies having star-formation rate between —0.5 < 
log[()(SFR/M., yr^') < 2. The grey-scale and contours are as for Fig.|5] and the range of logjQ(SFR/M3 yr"') and number of galaxies in each sub-sample 
are listed above each figure. 



We use p ri 

iBoroson & GreenI Il992h to quantify the relationships between 
these variables. This technique creates eigenvectors that are made 
up of linear combinations of the input variables, and are rotated 
to span the directions of maximum variance in the data. Principal 
component (PC) 1 contributes 71 per cent of the total variance, 
with weightings for the normalised SFR, mass, metallicity and 



extinction parameters as listed in Table[3]- i.e. the majority of the 
total variance in our sample can be explained simply by all four 
variables being correlated with each other, as has been found in 
previous sections. If all of the normalised quantities contributed 
equally to this PC, then they would all have weights of l/\/4 = 0.5; 
we find that mass contributes slightly more than metallicity or SFR 
to the first PC, but that all are broadly equivalent. 
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Figure 7. The variation in Ha extinction with star-formation rate and metallicity, for sub-.samples of galaxies having stellar mass between 8.5 < 
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Table 3. The four principal components, along with the percentage of vari- 
ance in the data that they explain. Parameters are as described in Table|2] 



0.5 





SFR' 


'Mass' 


'Metallicity' 


'Extinction' 


PCI (71%) 


0.469 


0.559 


0.453 


0.512 


PC2 (17%) 


-0.679 


0.011 


0.733 


-0.040 


PC3 (9%) 


0.313 


0.331 


0.237 


-0.858 


PC4 (3%) 


-0.471 


0.760 


-0.448 


-0.002 



PC2 contributes a further 17 per cent of the total variance, 
with significant weightings of —0.679 for SFR and 0.733 for metal- 
licity (see Table Isj - i.e. the majority of the remaining scatter in 
our data can be attributed to an anti-correlation between SFR and 
metallicity, for little variation in mass and extinction. A weak anti- 
correlation can be seen in each of the panels of Fig.|7](as illustrated 
by the contoured regions). However, we have not fully removed the 
effects of stellar mass in these panels, as there is still a 0.5 dex varia- 
tion in stellar mass across each of the sub-samples. Note that since 
the weightings of mass and extinction in PC2 are very small, the 
anti-correlation between SFR and metallicity will partially cancel 
out the correlation between SFR and metallicity from PC 1 . 

The direction of PC2 implies that, if the stellar mass depen- 
dencies of other variables are removed from our dataset, we should 
see an anti-correlation between SFR and metallicity. In order to re- 
move these dependencies, we bin the full dataset by stellar mass 
and fit a polynomial to the median mass and SFR for each bin, 
which gives 



logio(SFR/M0yr"') = O. 64 + 0.76 ^' + 0.12^2 



(3) 



where X = logio(M*/lOi"M0). 

We repeat this with the binned median mass and metallicity, 
to obtain 






'o 



-0.5 



""' " — -._^«i 




SFR residuals 



12 + logio(0/H) = 8.68 + 0.23 X-0.13X2 



(4) 



Figure 8. The distribution of SFR and metallicity residuals, after stel- 
lar mass dependencies are removed using equations (5) and (4). SFR 
residuals are given by logio(SFR/SFR„,o;]el). ^nd metallicity residuals by 
[12 + log[o(0/H)] - [12 + log|o(0/H)„io^ei]. The distribution of residuals 
is indicated by the grey-scale, and the solid and dashed lines indicate how 
the metallicity residual varies with SFR residual (median, ± Iff). 



In both cases, moving to a higher-order polynomial did not signifi- 
cantly improve the fit. 

We use these relations to predict the SFR and metallicity of 
all galaxies in the sample (the redshift range studied is sufficiently 
small that any cosmic evolution of these relations is negligible). By 
subtracting these predicted values from the observed values, we ob- 
tain the amount of SFR or metallicity variation that is not related to 
the overall stellar mass dependency, which we denote as the SFR 
and metallicity residuals. The relationship between SFR and metal- 
licity residuals is shown in Fig. [8] An anti-correlation between SFR 
and metallicity, after stellar mass effects are removed, appears in 
the data for positive SFR residuals. Note that PCA only finds linear 
relationships between variables - as we have had to fit quadratics 
in order to represent the data successfully, we do not recover a pure 
anti- correlation between the two residuals. 

Ellison et al.l ( 120081) have found a similar trend that for a given 
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stellar mass, galaxies with a high SFR have systematically lower 
metallicities than their low-SFR counterparts. They also find that 
more concentrated galaxies, with a smaller half-light radius, tend 
to be more metal-rich. Their interpretation of these results is that 
galaxies with a high surface mass density may have had a higher 
star formation efficiency in the past than their lower-density coun- 
terparts - this would lead to an increased production of metals, but 
lower current SFRs as more of the gas reservoir would have been 
depleted. This interpretation is consistent with our results. 

PCS contributes 9 per cent of the total variance, with the only 
significant weighting being for extinction. This may imply that 
there is a residual scatter in the extinction of each galaxy that is 
caused by some property other than mass, SFR or metallicity. PC4 
contributes the remaining 3 per cent of the variance, and will be 
dominated by noise in the sample. 



3.4 Modelling the variation in extinction 

Having shown that extinction is correlated with various properties 
of a galaxy, we now attempt to model these variations, and identify 
the amount of information that is required in order to accurately 
estimate the typical extinction of a galaxy. We will quantify the 
accuracy of our models in two ways: 

(i) calculating the residuals between the model and observed ex- 
tinctions for all galaxies in our sample, and measuring the width of 
the Gaussian distribution that describes these residuals; 

(ii) plotting the distribution of residuals against galaxy mass, 
SFR and metallicity, in order to identify any significant trends that 
have not been accounted for in the models. 

Our first model is the assumption that all galaxies can be rep- 
resented by a single value of extinction, which we take to be the 
median value of Ana in our sample, 1.03 mag. Fig. |9ta) shows 
that the distribution of residuals is roughly Gaussian, with width 
CT = 0.46 mag. The strong variations in extinction with SFR, metal- 
licity or stellar mass are shown in Figs. |9fb) to (d), as previously 
identified in Section [3711 We show both the uncertainty in the me- 
dian estimate, and the ±lcT distribution widths, in order to fully 
quantify the distribution. This model is the one most commonly 
used in the literature for Hot-selected samples, but is clearly not a 
good fit to the data. 

Our next three models come from polynomial fits to the binned 
data from Section lJH in the form of 



Table 4. Co-efficient.s of the polynomial fits to the extinction as determined 
using Eq.|5] 



A\ia ■ 



!=0 



(5) 



The extinction / SFR relationship, with X = logjo (SFR/Mq yr^ ' ), 
is well-fitted by a linear relationship. The extinction / metallicity 
relationship, with X = [\2 + log]o(0/H)] — 8.5, requires a much 
higher-order polynomial (5th order) in order to model the curva- 
ture at high metallicities. We note here that it is unlikely that an 
extinction / metallicity relationship would be of much practical use 
- if sufficient spectroscopy exists that a reliable metallicity can 
be calculated, it would be better to measure extinction values di- 
rectly. However, we include this model in our discussion for com- 
pleteness. Finally, the extinction / stellar mass relationship, with 
X = logio(M*/10 MfT,), can be modelled by a 4th-order polyno- 
mial. These relationships are overlaid on Fig. |4] for reference, and 
all are a good fit to the binned data - we do not obtain a significantly 
better fit by moving to higher-order polynomials. The co-efficients 
of the polynomial fits are given in Table ID 



Model 



Bo 



B2 



Bs 



54 



SFR 


0.53 


0.54 






Metallicit 


y 0.64 


1.75 


2.96 


-2.57 


Mass 


0.91 


0.11 


0.11 


-0.09 


ii 


1.11 


0.69 


-0.19 


-0.18 



-12.4 



The overall distributions of residuals are shown in Fig. |9le), 
(i) and (m), and again are consistent with being Gaussian, with zero 
mean and widths of 0.33, 0.33 and 0.28 mag for SFR, metallicity 
and mass. All three models show significant improvement over the 
assumption of a single extinction for all galaxies. The contribution 
to these residuals from uncertainties in measuring the Ha and H/3 
fluxes is about 0.2 mag, from Fig.|2lf), which is about half of the to- 
tal if errors are independent, and added in quadrature. Measurement 
errors in the SFR, metallicity and stellar mass parameters will ac- 
count for some of the remaining scatter, but are unlikely to account 
for all. This suggests that there remains up to 0.2 mag of intrinsic 
scatter between different galaxies (cf. PC3 in Section 3.3). 

The distribution of residuals for each of these models is shown 
in Fig.|9]as a function of SFR, metallicity and stellar mass, in or- 
der to test for remaining trends. By construction, the SFR residu- 
als show no trend when the SFR model is used, and equivalent re- 
sults are found for the other two models. However, the SFR inodel 
leads to a considerable variation in the residual extinction with ei- 
ther metallicity or stellar mass, in the sense that the model over- 
predicts the extinction of low-metallicity or low-mass galaxies. A 
similar result is found for the metallicity model, with extinction be- 
ing over-predicted for low-SFR or low-mass galaxies. However, the 
mass model shows significantly smaller residual trends, as expected 
from the results of Section. 3.2 and Figure 7. This implies that it de- 
scribes the data with less bias than the alternatives, although some 
weak trends with SFR and metallicity do still remain. 

We attempt to improve upon the mass model, by adding a 
further SFR-dependent term to equation ([Sj with coefficients of 
Bo = -0.09 and fii =0.13 (the best-fit values after taking out the 
mass dependence); this will produce an improvement if the scatter 
around the SFR-mass relation correlates strongly with extinction. 
The results are shown in the bottom row of Fig. |9]- a slight im- 
provement in the variation of residuals with SFR is seen (by con- 
struction), but the residuals show a worse trend with stellar mass 
than the previous model (since the strong SFR-mass correlation 
means that the additional SFR fit worsen the mass fit); the overall 
Gaussian width remains unchanged, at 0.28 mag. The largest resid- 
uals are seen at high stellar masses, possibly because this is the 
region where the mass / SFR correlation begins to break down, due 
to the effects of downsizing (where the most massive galaxies have 
already formed their stars, and undergo smaller amounts of current 
star formation than their less-massive counterparts). Concurrent fit- 
ting of the SFR and mass polynomials (even with the inclusion of 
cross-terms) was also attempted but does not decrease the overall 
Gaussian width of the residuals: although the trends of the resid- 
uals with SFR and mass improve slightly, those with metallicity 
worsen. Equivalent results are found if either a metal-dependent 
term is added to equation l|5j, or both a SFR-dependent and metal- 
dependent term are added, or all three terms are fitted concurrently 
- in all cases, the residuals are consistent with a Gaussian having a 
width of 0.28 mag. These results confirm that knowledge of stellar 
mass is sufficient to model the extinction of a galaxy in a statistical 
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Figure 9. The distribution of residual extinctions for our sample, for each of the models described in Section IT?) Each row represents a different model, which 
is named in the left-hand panel. The first column shows the overall distribution of residuals, along with the best-fitting Gaussian centered on zero. Column 2 
shows the variation in residuals with SFR, column 3 shows the variation with metallicity and column 4 shows the variation with stellar mass. Points with error 
bars indicate the median and error in the median; dashed lines indicate the ± 1 tJ distribution width. 



sense ~ once the mass variation is removed, very little improvement 
is obtained by a more complex model. 



Finally, we note that very recentlv lMannucci et alJfeOlOf) have 
proposed a 'fundamental metallicity relation', whereby galaxies 
define a tight surface in the 3-dimensional space of mass, SFR 
and metallicity. They introduce a new quantity, /x = log[()(Mass) 
- 0.32 logio(SFR), against which they argue that metallicity scatter 
of nearby galaxies is minimised. To test whether this can improve 
extinction estimates, we model the median extinction as a function 
of jl using a 4th-order polynomial (with co-efficients given in Ta- 
ble |4]l. However, we find that the resultant Gaussian width of the 
residuals relative to this model is 0.31 mag, with residual trends 
remaining with each of mass, SFR and metallicity that are worse 
than those of the mass model. We conclude that, for modelling ex- 



tinction, the stellar mass alone provides a better (and simpler) input 
than /i. 

3.5 Ha S/N bias 

One of our primary selection requirements was that sources have an 
Ha S/N ratio greater than 20, as measured from the observed line 
flux. By doing this, we are implicitly biasing our sample against 
highly extinguished sources - the error on the line flux measure- 
ment will be independent of the extinction of a galaxy, while the 
observed flux will be lower for galaxies with higher extinction. 

In order to confirm that this bias does not significantly in- 
fluence our conclusions, we repeat our analysis using an alterna- 
tive selection procedure which is less prone to an extinction bias. 
For galaxies in which the Ha and H/3 emission lines are both de- 
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tected with S/N ratio greater tlian 3, we use equation l|T} to calculate 
the Ha extinction. Similarly, for those galaxies with Ha S/N ratio 
above 3, but Hj3 S/N ratio below 3, we estimate the extinction us- 
ing the 3C7 limit on the Hj8 flux (strictly this is a lower limit on 
the extinction). We then correct the Ha line fluxes for the calcu- 
lated extinction, and derive the S/N ratio with which the Ha line 
would have been detected in the absence of extinction. We then set 
a primary selection requirement that this unextinguished Ha S/N 
ratio should be greater than 75, and we repeat the analysis of the 
paper. This selection procedure does not completely remove the ef- 
fects of extinction bias (the most extinguished sources may still be 
excluded if their observed Ha S/N ratio is below 3, or if their H/3 
line is undetected leading to an underestimated extinction and un- 
dercorrection of their Ha flux), but it does greatly reduce the bias. 
Therefore, if the same trends are found, that would provide confi- 
dence in their reliability. 

The median extinction of Ha in the new sample is 1.35 magni- 
tudes, significantly higher than the previous selection. This increase 
is driven both by the exclusion of low-extinction galaxies (mostly 
at lower masses; these were in the original sample but fail to reach 
the higher threshold for the revised sample) and by the inclusion of 
additional high-extinction galaxies (especially at higher masses). 
We re-calculate the relations with extinction determined in Sec- 
tion |3.4| The fitted polynomial parameters for the mass / extinction 
relationship areSo = 109, Si = 0.65, S2 = 0.12 andfis = -0.031; 
this has almost exactly the same shape as the previous relation, but 
is shifted up by about 0.18 magnitudes of extinction. Equivalent 
shifts are seen in the extinction relationships with other parame- 
ters. The distributions of extinction residuals (cf. Fig.[9) again show 
similar trends, with a scatter of 0.40 around the median extinction 
value being reduced to 0.37, 0.34 and 0.28 when taking out the 
SFR, metallicty and mass models, respectively. This confirms that 
mass is the predominant factor influencing the extinction of galax- 
ies; once again, adding in additional SFR- or metallicity-dependent 
terms does not further decrease the scatter. 

These tests leads us to conclude that although we are exclud- 
ing some of the more highly extinguished sources from our sample, 
the trends determined from the analysis are robust. 



4 DISCUSSION AND CONCLUSIONS 

We have used the superb spectroscopic coverage of SDSS DR7 to 
disentangle the effects of three of the physical parameters that are 
known to correlate with the dust extinction of a galaxy. We select a 
clean sample of star-forming galaxies using the BPT diagram, and 
calculate extinctions from observed Balmer decrements. We calcu- 
late metallicities from the 03N2 indicator, and after correcting for 
dust extinction and aperture size, we calculate SFRs from the Ha 
line. Stellar mass estimates have been calculated independently by 
the MPA/JHU group. 

We confirm that the typical dust extinction of a galaxy in- 
creases with greater SFR, metallicity or stellar mass. By selecting 
various sub-samples of galaxies, we show that the dominant factor 
is stellar mass - we conclude that the observed correlations be- 
tween extinction and SFR, or extinction and metallicity, are largely 
secondary effects brought about by the known dependence of these 
parameters on the stellar mass of a galaxy. The implications that we 
draw from our results are: 

• As galaxies build up their stellar mass over time, they also 
build up their dust content. This increase in dust leads to a greater 



Ha extinction, and therefore a dependence of extinction on stellar 
mass is seen. 

• Massive galaxies are able to retain a greater fraction of the 
metals that they produce, due to their deep potential wells, while 
less-massive galaxies lose metals from their shal lower potential 
wells through the effects of galactic winds (e.g. iTremonti et al.l 
l2004r) . A relationship between metallicity and extinction is there- 
fore formed - this is a secondary effect brought about by the de- 
pendence of both extinction and metallicity on stellar mass. 

• More massive galaxies are capable of forming a greater 
amount of stars in a starburst, and therefore a relationship be- 
tween the SFR and extinction of a galaxy is set up - this is also 
a secondary effect brought about by the relationship between SFR 
and stellar mass. This relationship will only hold for galaxies that 
are currently forming stars, and will break down at the highest 
masses, since these galaxies formed their stars at earlier epochs 
JCowie et alj|l99 6). However, extremely high-mass galaxies are 
poorly represented in our sample, since it was selected by the pres- 
ence of Ha emission to be undergoing current star formation. 

• An anti-correlation is seen between SFR and metallicity, for 
galaxies of a given stellar mass. This may be due to variations 
in star formation e fficiency between galaxies, as discussed by 
lEUison et alj OOOSh . where galaxies that have formed more stars 
in the past have a lower current SFR, and higher metallicity. 



Using equation ^ it is possible to predict the extinction of 
a galaxy with a given stellar mass, to within a typical error of 
0.28 mag (note that as discussed in Section 3.5, samples selected 
with different selection effects may require small offsets in nor- 
malisation of this equation, but the functional form and the scat- 
ter around it is largely robust to changes in selection technique). 
This typical error is broadly comparable to the accuracy with which 
extinctions can be estimated from the Balmer decrement, for our 
sample of galaxies selected with an Ha S/N ratio > 20, although 
this comparison will of course depend upon the S/N ratio of avail- 
able Ha and H/3 observations. We believe that the stellar mass of 
a galaxy is the best indicator to use when estimating the typical 
amount of dust extinction, both because it appears to be the funda- 
mental parameter that governs dust extinction, and also because it 
can be estimated in a way that is almost una ffected by dust, using 
the re st-frame A'-band luminosity alone (e.g. lLonghetti & Saraccd 
l2009r) . For some applications, it might be simpler to use the SFR- 
dependent correction (for example, correcting an observed Ha or 
UV luminosity function for the effects of dust). However, we cau- 
tion that this form of correction has the potential for introducing a 
significant mass-dependent bias into the results, and should be used 
with care. 

Estimating the extinction of truly 'passive' galaxies (with SFR 
< 0.1 M,;, yr^') is outside the scope of this work, as these are not 
detected in our Ha selected samples. However, for massive galax- 
ies, we note that the Ha selection limit corresponds to star forma- 
tion rates well over an order of magnitude below their historical 
average. These galaxies display high levels of extinction, implying 
that high-mass 'passive' galaxies are also likely to have high lev- 
els of extinction of the small amount of Ha emission that they do 
produce. 

We note that although all of our results have been derived 
for Ha extinction, they can be generally applied to line emission 
at another w a veleng th through the use of equation lUll and the 
ICalzetti et alj 1 I2OO O') dust attenuation law. If the extinction of UV 
radiation is required, rather than line emission, then a factor of 0.5 
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needs to be applied to this equation - see lCalzetti et alj ( 120000 and 
iGam et alj ( l2009bh for further details. 

Estimates of star-formation at high redshift (z > 1) are com- 
monly made using a wide range of different estimators: UV contin- 
uum; [Oil] or Ha emission lines; sub-mm, mid-IR or far-IR con- 
tinuum; radio emission. Each of these has their own advantages and 
disadvantages, with different sensitivities, angular resol utions, re- 
liabil ity of calibration, and effects of dust extinction (cf. lKennicutll 
1 1998!) . The latter effect is particularly important for star-formation 
rates estimated from Ha or UV selected samples (the most sensi- 
tive selection methods currently available); previous star-formation 
rate estimates from such samples have either required assumptions 
about the global extinction of all galaxies (e.g. a constant - a hy- 
pothesis which can be ruled out by the results in this work), spec- 
troscopic observations of each galaxy to measure individual dust 
extinctions (time-consuming, and technically challenging at high 
redshift), or a sufficiently large number of photometric data-points 
to enable SED-fitting techniques to estimate extinctions without 
issues of degenerate solutions. The results presented in this work 
suggest that a representative extinction correction can be derived 
from the mass, which can be estimated using far fewer photometric 
data points. Provided that the relation between mass and extinction 
is confirmed to be independent of redshift, this result could have 
important applications for studying large samples of star-forming 
galaxies, especially at high redshifts. 
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